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Abstract 

This paper deals with reduction of non-homogeneous linear systems of first order operator equations with constant 
coefficients. An equivalent reduced system, consisting of higher order linear operator equations having only one 
variable and first order linear operator equations in two variables, is obtained by using the rational canonical form. 
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1. Introduction 



Let K he a. field, V a vector space over K and \ei A -.V — > V he a, linear operator. Linear system of 
first order ^-operator equations with constant coefficients in unknowns xi, is: 



^(a;i) = 6iixi + 6122:2 + ■ . • + hinXn + ^\ 

A(x2)^h2\Xx + h22X2 + ■ . • + h2n,Xn + ^2 
^ A{Xn)^hn\Xx + hn2Xl + ■ ■ • + ferm^Jn + ^ 



(1) 



nxl 



for hij G K and tpi S V . We say that B = G iiT"^" is the system matrix, and <^= [<pi ... G V 

is the free column. If (^=0 system ([T]) is called homogeneous. Otherwise it is called non- homogeneous. 

Let x—\x\ ... Xr^ he a column of unknowns and A : — y ynxi j^g vector operator defined 
componentwise A{x) = [A{xi) . . . A{xn)]'^ ■ Then system ([1]) can be written in the following vector form: 

A{x) Bx + 0. (2) 

A solution of vector equation © is any column veV^^^^ which satisfies 
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Powers of operator A are defined as usual A* = A'^^o A assuming that A''^ = I : V — > V is the identity 
operator. By n*'* order linear A-operator equation with constant coefficients, in unknown x, we 

mean: 



A^{x) + biA''-\x) + ... + bnIix) 



(3) 



where bi, . . . ,bn&K are coefficients and ip^V. If (p = equation is called homogeneous. Otherwise it is 
non- homogeneous. A solution of equation ^ is any vector v^V which satisfies (jS)). 

2. Notation 

Let MgK^^" be an n-square matrix. We denote by: 
i^k) —^kiM) the sum of its principal minors of order fc (1 < fc < n), 

(^fc) '^fe =^fc(-^) the sum of its principal minors of order k containing i*^ column (1 < i,k < n). 

Let vi, . . . ,Vn be elements of K . We write M*(wi, . . . , t;„) G for the matrix obtained by substituting 

v=[vi ... Vn]^ in place of i*'* column of Af. Furthermore, it is convenient to use: 

(^KM; v) = 51{M- «i, . . . = 5l(M\v^,. . .,«„)). 
The characteristic polynomial As (A) of the matrix B^K'^''^ has the form 

Ab(A) = det(A/ - B) = A" + diA"-i + . . . + d„_iA + d„, 
where 4 = {-l)^5k{B), 1 < fc < n; see [2, p. 78]. 

The companion matrix of a monic polynomial p{X) = A" + di A"^^ + . . . + dn-iA + (i„ is the matrix 



Cn 





— (^n ~ — 1 








-d2 



The matrix C is in the rational canonical form if it is block diagonal, C = diag (Cp^ , Cp^ , . . . ,Cp^), 
where Cp. is the companion matrix of a monic polynomial Pi{X) = A"' + d^^iA"'"^ + . . . + cii_„._iA + di^m of 
degree ni — deg{pi) > 1, and pi \ Pi+i for all i such that l<i<k. 



3. Main Results 



The reduction of system ([2]) to a partially reduced system will be divided into two steps. 
Firstly, by using some basic properties of the rational canonical form, the initial system may be transformed 
into a proper system for further study. 

Theorem 3.1 [On rational canonical form [1, p. 397]] Any square matrix B gK"^"^ is similar to a unique 
matrix C in the rational canonical form. The matrix C is called the rational canonical form of the matrix B. 

Lemma 3.2 Let C he the rational canonical form of the matrix B , i.e. there exists a regular matrix P such 
that C = P^^BP. Then the system given in the form A{x) — Bx + (/?, can he reduced to the system 



A{y)=Cy + r^, 

where ijj = P^^ip is its free column and y = P~^x is a column of the unknowns. 



(4) 
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The following lemma will be very useful in the second step of our reduction process. 



Lemma 3.3 // the matrix M has the form: 



M = 



bi 1 
62 



5„_i ... 1 
bn Cln-1 a„_2 . . . ai 



then it follows 



fc-i 



6l{M) = (-!)'=( ^6,a,-,-fefe) (1 < fc < n). 



Proof. Observe that the function is linear with respect to the first column of its matrix argument. The 
first column [61 ... &„]'^ of the matrix M is equal to bjej, where ej denotes j*'' column of the identity 

matrix of order n. Therefore S].{M) = J2^=i bj^kiM; ej). It remains to prove that 



5l{M-e,) = < 



0, k< j; 

-(-1)', fc-j; 

{-l)''ak-j, k>j. 



Let ttfej be an arbitrary nonzero minor in the sum dk{M;ej). Since it is a nonzero principal minor which 
contains 1 from the first column, it must include 1 from the first row. Hence, minor akj must also include 
the element 1 from 2"'' row and 3'''^ column. In the same manner we conclude that our minor contains 1 
from z*'' row and + column for i G {2, . . . ,j — 1}. Accordingly, minor akj must contain the first j 
columns and rows. 

For k<j there is no such nonzero principal minor. 

If k — j, then the only nonzero minor in the sum dk{M; ej) is the minor formed from the first j columns and 
rows. Expanding it by the first column, we easily conclude that its value is (—1)'^^^. 

For k> j the same reasoning implies that ak.j must contain all columns and rows from 1"* to j"*. Now, 
let / + 1 be the smallest index greater than j such that ak.j contains (Z + 1)*^ column. Consequently, the 
minor must include element a„_i. We apply zig-zag argument again to deduce that our minor has a block 
formed from the last n — l columns and rows. So, the only minor having the desired properties is the one 
obtained by deleting all columns and rows between (j + l)*'' and Z*''. It follows that k = j + n — I and 
ak,j = i-l)'-^ ■ (-l)"-'-ia„_z = (-l)'=afc_,. □ 

Remark 3.4 The matrix M is the transpose of the doubly companion matrix introduced in [4], see 
also [5-7]. 

We will now prove the main theorem dealing with a special case when the rational canonical form of 
the system matrix has only one block. 

Theorem 3.5 Assume that the rational canonical form of the matrix B is exactly the companion matrix C 
of the characteristic polynomial As (A) = A" + diA"^^ + . . . + (i„_iA + dn — Ac (A). Then linear system of 
first order A-operator equations ([T|) can be reduced to the partially reduced system: 
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fe=i 



(Tec) <^ 



2/2 
2/3 

2/« 

where the columns y = [yi 



A{y2) - -02 

. . . y^Y <^'>T-d V' = [i'l 



(5) 



ipn]'^ are determined by y — P i ip = P ^ip for 
a regular matrix P such that C — P^^BP. 

Proof. Applying Lemma [3^ system ([2]) is equivalent to system (|4]), which can be presented in the following 
form: 

A{yi) = 2/2 + ^^1 

^(2/2) 2/3 + 0^2 

Myn-l) = yn + V'n-l 

A{yn) ^-dnyi-dn-iy2- ■■■-diyn + ipn ^ 



> . 



(6) 



(7) 



We now consider the linear operator: 

Ac(A) = A" + + . . . + dn-iA + drJ. 

By eliminating ?;2, ■ • ■ , 2/n from the previous system, we obtain the operator equation: 

Ac(A)(yi)-(^""'(^i)) 

+ [A''-\ij2) + d,A^-\^j^)) 
+ (^"-3(V^3) + diA"-3(V;2) + d2A"-3(^^)) 

n 

+ dlV'n-l + d2'lpn~2 + ■ ■ ■ + dn-2lp2 + ^n-lV'l ■ 

Using formulae ([S]), © and Lemma we complete the proof. □ 

Remark 3.6 It is important to say that the reduction from the previous theorem is reversible, i.e. reduced 
system ([5]) implies system ([T]). 

Remark 3.7 System ^ is called partially reduced because it has only one variable separated from the 
others by a linear n^^ order A-operator equation. 

We now continue with the general case of partial reduction. 

Theorem 3.8 Assume that the rational canonical form of the matrix B is the block diagonal matrix 

C — diag(Ci, C2, . . . , Ck) for some k, 2 < k < n, 

where Ci are the companion matrices of the monic polynomials Ad (A) = A"' +dj^iA"'~^ + . . .+c?i.„._iA+di.„. 
of degree Ui > 1 and Ac. | A^^i for all i, 1 < i < k. 

Let £1 = and let £i = i, 2 <i < k. Then linear system of first order A-operator equations ([T]) 

can be transformed into the conjunction of the partially reduced systems: 

k 

A (T^c.), (8) 



where every subsystem {TZd ) is of the form: 

rii 



yu+2 
y£i+3 

^ and ijj -■ 



k=l 

A{yii+i) - V'^i+i 

A{yii+2) - tpei+2 

[V'l . . . ipn]'^ o,i^e determined by y = P^^x and ip — P^^ip for a regular 



Columns y = [yi ... j/; 
matrix P such that C~P^^BP. 

Proof. According to Lemma [3.21 system ([2]) is equivalent to the following system: 

/\ A{y^) = ay, + ip. 



(9) 



where yi = [yi^+i 



yu+nA 



is a column of unknowns and V'i = H'ii+i ■ ■ ■ i^ii+ni]'^ is the free column of 



the linear system A(yi) = Ciiji + V'i. We now consider the linear operator: 

Ac.(A) = A"' + + . . . + d,,„,_iA + d,,„,/. 

Therefore, by Theorem 13.51 every subsystem A{yi) — Ciiji + is equivalent to the partially reduced system 



{TtCi ) . We now obtain the equivalence between system ([9|) and the system /\ [TZd ) ■ □ 

i=l 

k 

Remark 3.9 System /\ {TZci) is equivalent to the initial system ([1]). Every subsystem {TZci) consists of a 

1=1 

higher order linear A-opcrator equation, in only one unknown, and first order linear A-operator equations, 
in exactly two variables. 

Remark 3.10 The initial system has a solution if and only if every higher order linear A-operator equation 
of the partially reduced system has a solution. 

Remark 3.11 The above mentioned reduction formulae can be applied to two special types of non- 
homogeneous systems, to systems of differential and systems of difference equations with constant coef- 
ficients. The usual way of proving the reduction formulae for a non-homogeneous system of differential 
equations with constant coefficients is to use the Jordan canonical form. However, this method is somewhat 
complicated, see [3, pp. 166-178]. We give a better and more efficient alternative formulae of reduction for 
an arbitrary A-operator. Our method involves the rational canonical instead of Jordan form which makes 
reduction process much more easier to handle and enables us to find general solutions. 

4. Examples 

The reduction formulae in special cases of n = 2 and n = 3 will be given as examples. 
Example 4.1 Consider the linear system of the first order A-operator equations, in unknowns Xi,X2- 



A{xi) = 6iia;i -I- bi2X2 + (fii 

A{X2) = b2lXi + b22X2 + "^2 



(10) 



Depending on the block structure of the rational canonical form C of the matrix B, partially reduced system 
of system (ITUl) has two possible forms. 
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1° The rational canonical form C of the matrix B has one block. Then the partially reduced system has 
the form: / , , 

As(A)(yi) - 5liC;Ai^,),Ai^P2))-5liC;^Pu^2) 

y2 = A{yi) - -01 

2" The rational canonical form C of the matrix B is diagonal matrix C — diag(Ci, Ci). Then system ([TO|l 
is equivalent to the partially reduced system: 

Aci(A)(yi) -01 
Aci(A)(y2) -^2 

Example 4.2 We now deal with the system of first order A-operator equations in three unknowns xi,X2,x^: 

A{xi) = biixi + 612X2 + 6i3a;3 + <Pi I 

A{X2) ^ b2lXi + b22X2 + b23X3 + ^2 )■ (11) 
A{x'i) = 631X1 + 632X2 + 633X3 + (^3 ^ 

Partially reduced system ([8]) of system (ITlT) depends on the block structure of the rational canonical form C 
of the matrix B. The following cases are possible. 

1° The rational canonical form C of the matrix B has one block; then © has the form 

3 

AB(A)(yi) = Y.^-lt^^6\{C-A^-\i;^),A^~Hi^2),A^-\i:^)) 
fe=i 

y2 = A{yi) - -01 
^ ys = A{y2) - -02 

2" The rational canonical form C of the matrix B has two blocks C = diag(Ci, C2), where Ci and C2 
are matrices of order 1 and 2 respectively; then (|5]) has the form 

Ac,(^)(yi) = 0i ] 
^cM){y2) = 5\{C2;A{rP2),A{^^)) - 5^(C2;^2,^3) 
yz = A{y2) - 02 

3° The rational canonical form C of the matrix B has three blocks equal to Ci ; then ([8]) has the form 

Ac,(^)(yi) = 0i 1 
Aci(^)(y2) = 02 
Aci(^)(y3) = 0^3 
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